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Abstract 

The mechanism for the relaxation of the cosmological constant is studied and 
elaborated. In the model used for the analysis of the relaxation mechanism the 
universe contains two components: a cosmological constant of an arbitrary size 
and sign and a component with an inhomogeneous equation of state. Owing to the 
dynamics of the second component the universe asymptotically tends to a de Sitter 
phase of expansion characterized by a small effective positive cosmological constant. 
An analysis of the asymptotic expansion for a general inhomogeneous equation of 
state of the second component is made. Several concrete examples are presented 
and the stability and speed of convergence to their fixed points are analyzed. It is 
found that the speed of convergence to a fixed point is large whenever the absolute 
value of the cosmological constant is large. 



1 Introduction 



The understanding of the structure, dynamics and composition at cosmological 
scales is one of the biggest challenges for physics and maybe even for science in 
general. The results of cosmological observations in recent years [H El [3] and nu- 
merous attempts of their theoretical explanation [4J have lead us to an astonishing 
picture of the present universe. We presently understand only a minor part of the 
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composition of the universe whereas its unknown part is attributed to the cos- 
mic dark sector: dark matter and dark energy. Additional pecuhar feature of the 
universe at the present epoch is its accelerated expansion. The dynamics or mecha- 
nism leading to the accelerated expansion became the hot topic in cosmology soon 
after the accelerated expansion had been established observationally . Almost ev- 
ery element of our, until then valid, picture of the universe and its dynamics, came 
under scrutiny in search of the acceleration mechanism that would be consistent 
with the observational data. 

First attempts intervened into the composition of the universe: the existence 
of a new component with negative pressure (or even several of them) was postu- 
lated. The proposals for the nature of this new component, called dark energy, are 
numerous [5- ^o name just a few, dark energy models range from the dynamics 
of scalar field in the potential [5], which further comprises quintessence, k-essence 
and phantom energy models, various cosmic barotropic fluid models [S] to the very 
concept of static or dynamic cosmological term [7]. Further attempts towards the 
acceleration mechanism questioned general relativity as a theory of gravitation at 
cosmological scales p]. The very concept of our universe as four-dimensional was 
also relaxed in braneworld models in the quest for the explanation of acceleration 
[9]. The list of ingenious approaches to this problem certainly does not end here. 

The problem of the cosmological constant (CC) |10[ \TT\ [T2] enjoys a special 
place among the multitude of proposals for the acceleration of the universe. The 
ACDM model is probably the simplest and the most popular model of the dark 
sector. The CC problem is an old one and it comes in many forms. From the 
viewpoint of fundamental theories, such as quantum field theory (QFT), the size 
of A consistent with observations is in a notoriously huge discrepancy with the 
predictions of these theories. This is a so called "old CC problem". An additional 
problem comes from the fact that the present energy densities of matter and CC 
are of the same order of magnitude, despite their significantly different scaling 
with the expansion of the universe. This is a so called "coincidence problem". 
Although the ACDM model generally fits the data very well, other challenges to 
ACDM cosmology have been identified [13j- However, the importance of the CC 
problem and its possible resolution is even bigger given that it underlies many 
other models of the cosmic acceleration, such as dynamical dark energy models. 

In this paper we further elaborate the model of relaxation of the cosmological 
constant introduced in and studied in [15] . The relaxation of the cosmological 
constant is understood as a dynamical process in which the universe with a large CC 
asymptotically tends to a de Sitter regime. The Hubble constant characterizing 
this de Sitter phase is equivalent to a small effective cosmological constant, = 
Ae///3. We study the stability and convergence properties of the said model and 
provide several examples. 

2 The model of the cosmological constant re- 
laxation 

In this section we present a short summary of the CC relaxation model introduced 
in [T3] and discussed in [15]. We consider a two component cosmological model. 
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The first component is a cosmological constant of arbitrary size and sign, wliereas 
the second component is described by an inhomogeneous equation of state (EOS). 
The expansion of the universe is given by the Friedmann equation 

H^ = '-^{P^ + P). (1) 

The dynamics of the second component is described by the inhomogeneous EOS 
of the type 

p = wp- 3Co^"+' . (2) 

The formahsm of the inhomogeneous EOS proves to be very useful in the de- 
scription of various cosmological phenomena [161 IIZl ISl IHl [201 [211 [22] . It has 
been shown [IB] that the inhomogeneous equation of state can be interpreted as 
an effective description of modified gravity or braneworld dynamics |23l . An- 
other possible interpretation of an inhomogeneous EOS related to bulk viscosity 
[251 [261 [271 [281 [291 [301 [31] and in particular its generahzations. 

Using ([T]) and ^ with the standard continuity equation for the second compo- 
nent we obtain a dynamical equation for the Hubble function H^: 

dH' + 3{1 + w)^(h'- 5^ - ^(//2)(a+i)/2\ ^ . 

a \ 3 1 + w J 

Using the notation 



h={H/Hxf, s = a/ax, X = SnGpA/SHj,, ^ = SttGCqH^^^ /{I + w) , (4) 
where H{ax) = Hx, we arrive at the equation relating dimensionless quantities 

s— + 3(l + u;)(/i- A-C/i("+^)/^) = 0, (5) 
as 

and the initial condition /i(l) = 1. 

The relaxation mechanism for the cosmological constant is realized in the regime 
a < — 1. For an analytically tractable case a = —3, it is straightforward to show 
[14] that for large absolute values of the rescaled CC term A of both signs it is 
possible to obtain the relaxation of the cosmological constant. 

For negative values of A and A^ » ^ with > and w > —1, the asymptotic 
expansion of the universe is given by a small value of the Hubble function 

{l+w)\k\ {l + w)\pjy\ 3 ^ ' 

This small value of corresponds to a small effective positive cosmological con- 
stant A~jj. A crucial result is that the Ajy^ is small because A is large in absolute 
value, as presented in ([6|). The dynamics of the Hubble function for negative A 
and some typical values of parameters is given in Fig. [H 

For large positive values of A (so that A^ ^ |^|) and for ^ < and w < — 1, the 
Hubble function tends asymptotically to 
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Figure 1: The dynamics of the scaled Hubble function for parameter values a = —3, 
A = —5000, ^ = 0.03 and w = —0.9. The dynamics consists of two phases of accelerated 
expansion connected with an abrupt transition. 



,i„^^ = ^J, = _i±%,__J<2_ = ^. (T) 

Again as in the case of negative A, the small asymptotic value of the Hubble 
function can be interpreted as a small effective positive cosmological constant 
A+ ^. The effect ive CC is small because A is large. The dependence of the Hubble 
function on the scale factor for positive A and some typical values of parameters is 
depicted in Fig. [2j 

As described above, a simple model defined by ^ and ([2]) provides parameter 
regimes in which the relaxation of the cosmological constant can be realized. It is 
important to stress that in this model there is no fine-tuning and the asymptotic 
de Sitter phase is characterized by an effective positive CC which is small because 
the size of the real CC |A| is large. Based on the calculations from e.g. QFT we 
could say that the size of A is "naturally large". Additionally, since ^ describes 
the deviations from GR or the generalized effects of viscosity we could say that 
is "naturally small". Then, from the behavior of the model for a < —1 p3J and 
in particular from ([6]) and d?]) we can expect that the size of the effective CC in 
the asymptotic de Sitter phase is "naturally small". This fact corresponds to the 
resolution of the old CC problem in the studied model. 



3 Stability and convergence analysis 

Given the potential of the relaxation mechanism described in the preceding sections 
for the resolution of the cosmological constant problem, it is important to study 
more general inhomogeneous equations of state for which the relaxation mechanism 
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Figure 2: The dependence of the Hubble function on the scale factor for the values of 
parameters a = —3, A = 5000, ^ = —0.03 and w = —1.1. Two phases of accelerated 
expansion are interconnected with an abrupt transition. 

can be realized. Furthermore, as the dynamics of the Hubble function in the model 
|14j exhibits abrupt transition followed by a swift stabilization at the asymptotic 
value, additional insight into the dynamical details of the approach to the asymp- 
totic value of H is needed. To this end, we consider a general inhomogeneous 
equation of state 



p = wp- g{H'^ 



(8) 



which, using the standard rescaling as given in dH), results in the following dynam- 
ical equation for the Hubble parameter: 



dh 

s— + ?,{l + w){h 
as 



A-/W) = 0, 



(9) 



with the initial condition /i(l) = 1. Here we have /(/i) = {{8iTG)/{3{l+w)H]^))g{H]^h). 
We further introduce the notation F(h) = 3(1 -|- w){h — A — f{h)) and F'{h) = 
3(1 + _/'(/,)). 

For a relaxation mechanism to be effective, we expect the dynamical system 
([9]) to have a fixed point /i* which is much smaller than |A|, i.e. for which we have 
/i* -C |A|. From the condition F{h^) = and the said expectation of the size of /i* 
it is straightforward to obtain 



K = r\-^)- 



(10) 



The sign of the quantity F'{h^:) determines the stability of the fixed point /i*, 
whereas its size controls the speed of convergence to the fixed point. The fixed 
point /i* is stable for F'{h^) > and it is given by the expression 
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F'{K) = 3(1 + w){i- nr\-\)) = 3(1 + w)[i- (yzry^z^ j • (n) 

We further apply the general analysis displayed above to specific examples of 
the function f{h). First we turn to the choice f{h) = ^//i, already analyzed in 
detail in [TJ] and presented in the preceding sections of this paper. It is easy to 
see that in this case we have 

K = -| (12) 



and 



F'{h)=3{l + w)(l + ^] . (13) 



/i2 

Combining (fT2]l and (fT3]) and using ^ |^|, we obtain 



F'{h,) = 3{l + w)^. (14) 



This results sheds additional light on the results obtained in |14j and [15]. From 
()12p and (I14p it becomes clear that to have a stable fixed point for A > we must 
have ^ < and u; < — 1, whereas for A < we need ^ > and w > —1. However, 
the most important feature of (|14|) is that the speed of convergence towards the 
fixed point /i* depends quadratically on A. This finding explains why the fixed 
points in Figs [T] and [5] are reached so swiftly after the transition. 

In the following example we consider the function f{h) = Ai In h. For this 
functional form of f{h) we have 

K = (15) 

and 

F\K) = 3(1 + w)[l- ^le^/^i) . (16) 

The condition of the stability of the fixed point /i* is —3(1 + w)Ai > (given 
that the second term in (jl6p dominates). The speed of convergence to the fixed 
point grows exponentially with A. 

The next example employs the functional form f{h) = A2e^^^'^. As in the 
preceding examples we have 

K = (17) 



and 




F'iK) = 3il + w) l--{\n{-—)) . (18) 



The conditions for the existence of small and positive /i* comprise —X/A2 < and 
61/ ln(— A/yl2) > 0. The fixed point h^^ is stable for —(1 + w)X/bi < 0. 
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Finally, the last example to be considered in this paper is specified by the 
function f{h) = exp |^exp (x)) • further obtain 

= (19) 



In (in ( 



A, 



and 




F'(/i.)=3(l+«;) 1-^ln -- In In . (20) 



The expressions (|19p and (|2Up show the conditions for the existence of a positive 
fixed point K are < 0, ln(-A/A3) > and 62/ ln(ln(-A/^3)) > 0. The fixed 

point is stable if (1 + w)^ In [^-^ < 0. 

The conclusions of the general analysis given at the beginning of this section and 
the specific results for the studied examples indicate an interesting characteristic 
of the speed of convergence to the fixed point h^,. The faster the growth of the 
function f{h) when h acquires small values, the smaller the speed of convergence to 
the fixed point, measured by the dependence of F'{h^) on A. However, even for the 
examples with the mildest dependence of F'{h^) on A, the speed of convergence is 
extremely large owing to the fact that |A| is large. The examples discussed in this 
section were selected primarily to better illustrate this property of the convergence 
to the fixed point. This is especially true for the last example of double exponential. 



4 Conclusions 

The results of this paper, along with the findings of [H] and p^, show that the 
mechanism of the CC relaxation based on the inhomogeneous EOS is robust. In 
section [3] general conditions for the onset of the CC relaxation mechanism are given. 
Various inhomogeneous EOS can reproduce the asymptotic de Sitter phase without 
particular fine-tuning. A general needed feature of the inhomogeneous EOS is that 
the inhomogeneous term becomes increasingly important as the expansion slows 
down and its effect finally equilibrate the action of a large CC. Thus the expansion 
of the universe settles down in a de Sitter phase characterized by a small effective 
positive cosmological constant. An interesting feature of the models is that the 
faster the inhomogeneous part, defined by f{h) grows as h becomes small, the 
slower the approach to the asymptotic value, i.e. fixed point. However, for a large 
|A|, even in the examples with the slowest approach to the fixed point, the speed 
of convergence is large, principally because of the very size of |A|. These results 
largely explain some features of plots given in Figs [T] and [2j Namely, given a very 
large speed of convergence to the fixed points, it is easier to understand the abrupt 
transition in the dynamics of h and especially an extremely fast approach to the 
asymptotic value. All these results add to our understanding of the CC relaxation 
mechanism based on the inhomogeneous EOS and provide a further incentive to 
incorporate the said mechanism into a complete and workable cosmological model. 
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